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Abstract 



We study numerically the decay of massive and massless inflatons into massive 
excitations, via a (f) 2 X 2 coupling, in the expanding Universe. We find that a 
wide enough resonance can survive the Universe expansion, though account for 
the expansion is very important for determining precisely how wide it should 
be. For a massive inflaton, the effective production of particles with mass ten 
times that of the inflaton requires very large values of the resonance parameter 
q, q ^ 10 8 . For these large q, the maximal size of produced fluctuations is 
significantly suppressed by the back reaction, but at least within the Hartree 
approximation they are still not negligible. For the massless inflaton with 
a A(/> 4 /4 potential, the Universe expansion completely prevents a resonance 
production of particles with masses larger than \/A0(O) for q up to q = 10 6 . 
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The quantitative theory of reheating after inflation is an important problem in inflation- 
ary cosmology because the way the reheating proceeds determines possibility of certain 
phase transitions and scenarios of baryogenesis. It was recently realized that stimulated 
processes, due to large occupation numbers of produced bosons, can be most significant 
during reheating [0-[5|]. However, despite considerable progress made still little was 

known quantitatively about the important case when the coupling of the inflaton field (4>) 
to bosonic fields (X) is relatively large (the wide resonance case). It was not clear precisely 
how wide the resonance should be (for a given mass of X) to overcome the expansion of the 
Universe. If the coupling required by this condition turns out to be too large, the resonance 
will be prevented by non-linear effects, which in this case will set in when X fluctuations are 
still of the quantum strength. Thus even the existence of a wide resonance in the expanding 
Universe was not firmly established, let alone its quantitative details. The purpose of the 
present work was to study the wide resonance case by numerical means. 

The main quantity that one wants to know for a given inflationary model is the maximal 
size of fluctuations of fields produced in the inflaton decay. The size of these fluctuations 
determines which phase transitions took place after inflation |8|,|9j , whether enough of heavy 
leptoquarks could be produced to generate the baryon asymmetry [p|,|9|,[7l,|lQ| , and how strong 



was the cosmological breakdown of supersymmetry [pT| , p^| . 

The reason why fluctuations produced at reheating could be sufficiently large for the 
above effects to exist is that, as now recognized, reheating in many cases starts with a stage 
of parametric resonance (called preheating in ref. ||), during which the bosonic fluctuations 
grow extremely (exponentially) fast. In the expanding Universe, a parametric resonance 
is effective only when it is sufficiently wide or when the system is (classically) conformally 
invariant. Otherwise, the redshift of momenta prevents the resonance from developing [JT3j . 



The conformally invariant case is realized when a massless inflaton decays into a massless 
fluctuation field. Production of a massive field can only take place in the wide resonance 
regime. In that case, there can occur an interesting and somewhat counterintuitive phe- 
nomenon, when quanta of a field that strongly couples to the inflaton can be efficiently 
produced even if their mass is much larger than the frequency of the inflaton oscillations. 
The effect can be thought of as a "many— > few" process, in which many quanta of the 
inflaton field assemble together to produce a few heavy quanta. If the fields that can be 
efficiently produced by this mechanism include leptoquarks of a grand unified theory (GUT), 
the subsequent non-equilibrium decay of these leptoquarks could give rise to a sizable baryon 
asymmetry [gggj^. 

In the present work, we wanted, first of all, to find the regions of parameters, for which a 
wide resonance can overcome the redshift of fields and momenta in the expanding Universe. 
We considered two models: in the first one, the inflaton potential was m 2 2 /2, in the second, 
it was A0 4 /4. In both cases, we considered decay of the inflaton into a heavy scalar field 
X of mass M x , due to coupling ^g 2 (f> 2 X 2 . The width of the resonance in these models is 
determined by parameter q: in Model 1 (a massive inflaton) q = g 2 (p 2 (0)/Am 2 , where 0(0) is 
the value of the inflaton field when it starts to oscillate;[] in Model 2 (the massless inflaton), 
q = g 2 /AX. 



For a static universe, in the linear regime with respect to fluctuations, this agrees with the 
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In Model 2, we have found no resonance for any Mx > \/A0(O) with q up to q = 10 6 . 
In Model 1 with Mx/m = 10, which may be appropriate for scalar leptoquarks, we have 
found that the redshift is overcome only for very large q, q ^ 10 8 . Thus, a wide enough 
resonance can indeed survive in the expanding Universe, though account for the expansion 
is very important for determining precisely how wide it should be. 

Our second objective in this work was to determine the maximal size achieved by X 
fluctuations in Model 1. In general, both the Universe expansion and the back reaction 
play a role in determining the maximal size of fluctuations. A fully non-linear calculation 
along the lines of ref. which would allow us to include consistently all back reaction 
effects, was technically unfeasible for q 3> 10 3 . So, we took the back reaction into account 
approximately, using a simple Hartree procedure 

The Hartree approximation may overestimate the maximal size of fluctuations. The 
interaction term in the equation of motion for <fi becomes important when (X 2 ) ~ </> 2 (0)/g 
(the brackets denote averaging over space). One might expect that this value signifies the 
end of the resonance stage. However, in the Hartree approximation, when this value is 
reached, (X 2 ) merely renormalizes the mass of <fi 0] and continues to grow. We cannot 
say at present if that further growth is an artifact of the Hartree approximation in this 
model and will disappear in the fully non-linear picture. In the conformally invariant case 
(see below), the limited set of data points that we have for the full problem does fit on a 
1/q curve, in disagreement with the Hartree approximation. We are not willing to make 
any inferences for the present case, however, because the time evolution of is rather model 
dependent. 

We neglected any self-coupling of the field X. In the Hartree approximation, an estimate 
shows that a quartic self-coupling AX 4 /4 can be neglected during the entire resonance stage 
when A ^ g 2 . If the maximal size of fluctuations is in fact smaller than obtained in the 
Hartree approximation, this self-coupling can be neglected at the resonance stage even for 
larger A. 

The resulting (X 2 ) max as a function of q (the brackets denote averaging over space) is 
shown in Fig. |l|. It is suppressed for large q but is still not negligible. If we measure the 
strength of the (massive) fluctuations with regard to symmetry restoration |8]|| using the 
effective "temperature": T e 2 ff /12 = (X 2 ) max , then for q = 10 8 , the Hartree approximation 
gives T eff of order 10 _3 M P1 . This is of order of the GUT scale, so the process could be 
relevant for symmetry restoration even in Grand Unified models. 

One should keep in mind, though, that the size of fluctuations can be further altered 
(apart from the redshift) by the subsequent chaotic evolution ||, absent in the Hartree 
approximation. That chaotic evolution is rather fast because it is due to stimulated, i.e. 
classically enhanced, processes. For example, in applications to GUT baryogenesis, the 
chaotic evolution will in large part take place before the B-non-conserving decays of X par- 
ticles. Therefore, some quantitative understanding of it is required before reliable estimates 
of the baryon asymmetry can be made; this will be a subject of a forthcoming paper. 

Below we describe our models and calculations in more detail. 



definition of the parameter q of the canonical Mathieu equation [14 
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Model 1: Massive inflaton. In conformal variables, and after a convenient rescaling, the 
action for the matter fields is 
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where r = ma(0)r], r] is the conformal time; £ = ma(0)i; y? = 0a(r)/0(O)a(O), x = 
Xa(r)/0(O)a(O); a(r) is the scale factor of the Universe. 

We take r = at the moment of the first extremum of 4>{t). Solving the equation of 
motion for the zero-momentum mode of the inflaton coupled to gravity, together with the 
Einstein equation for the scale factor, we have found 0(0) = 0.28Mpi. At later times, we 
continued to determine the evolution of the scale factor self-consistently, including the in- 
fluence of produced fluctuations in the Einstein equations. After r = 0, a matter dominated 
stage quickly sets in, for which one can use 



(2) 



However, as our self-consistent solution confirms, for some values of the parameters the 
expansion later crosses over into a radiation dominated stage. That happens when the 
interaction potential \g 2 4> 2 X 2 overtakes the mass term of 0. 

In the Hartree approximation, the zero-momentum mode is the only non- vanishing mode 
of the field ip; we keep for it the notation 93. Its equation of motion is 

V ? --V ? + ^(^ + 4g(x 2 )^ = 0, (3) 
a Gr(0) 

where the brackets denote an average over space. This should be solved with the initial 
conditions ty?(0) = 1, 0(0) = 0. The equation of motion for the Fourier components of the 
field x is 

Xk + uJ 2 k (T) Xk = , (4) 

where 

We will sometimes use the notation 



M 



The Fourier components are normalized in such a way that the variance of the fluctuation 
field is 
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Because fluctuations are produced mostly with relatively small values of (comoving) mo- 
menta k, the inevitable presence of a momentum cutoff in our calculation, somewhere at 
larger momenta, will automatically result in a properly renormalized (x 2 )- If we wanted 
to embed our calculation into a fully quantum-field-theoretical picture, we would identify 
the parameters we use (m, g 2 etc.) with the running parameters of that picture, taken 
at a normalization point of the order of our momentum cutoff; moderate changes in the 
normalization point will only cause corrections of relative order g 2 /4:7t = qm 2 /4tt<J) 2 (0). 

The initial conditions for XV are randomly distributed with the probability density 



^fc[Xk] K ex P 



20 2 (O) 
m 2 



"*(0)|xk(0)l 



(8) 



According to eq. (|1|), m 2 /0 2 (O) is the parameter that regulates the strength of quantum 
fluctuations compared to classical ones in this model. The initial "velocities" are locked to 
their corresponding "coordinates" 

X k (0) = H^(0) + M0)] Xk , (9) 

where h = a /a. For a derivation of these initial conditions, see ref. f|. The semiclassical 
description is reliable as long as g 2 /4ir = qm 2 /4ir<j) 2 (0) <C 1. 

The results of numerical calculations for the realistic case m = lCT 6 Mpi are presented in 
Figs. 0-1 

Three factors can potentially suppress a parametric resonance in the expanding Universe, 
each represented by the corresponding term in the equations of motion. In our conformal 
variables, they are: the time dependent mass term oc M\ in eq. (|5]), which redshifts the 
modes away from the resonance; the time dependent mass term of the field tp in eq. @, 
which makes the oscillations of the background field aperiodic; and the back reaction term 
oc q in eq. (|3|). A very useful quantity to introduce is the effective q, 

where is the (decreasing) amplitude of the physical field 0, and m e s is that field's effective 
mass including the Hartree correction. In a slowly expanding Universe (the adiabatic limit), 
Qes{T~) would indeed be the correct resonance parameter at time r. 

To understand the effects of the redshift and the back reaction, it is convenient to see 
first what happens when the back reaction is neglected. This is indeed the correct limit for 
sufficiently small q (a more precise criterion depends on m x ). Then, m e s = m and 

For a parametric resonance to be effective at time r, q^{r) should not be much smaller than 
unity. Because a resonance will need some time tq to start, the condition for it to begin at 
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all, g e ff(7o) ~ 1, already requires large values of q = g e fr(0). By our definition, a resonance 
starts when (X 2 ) rises above quantum fluctuations. E.g. in Fig. [], tq ~ 6; some fc-modes 
can start to grow earlier, however. 

Consider first the case of massless X, Mx = 0. We have found that even then the 
resonance starts only for q ^ 10 3 , see Fig. [I]. With q ^ 10 3 , the parametric resonance 
starts at some point but terminates when q e g drops below unity. Hence, in the expanding 
Universe, the resonance stops and the maximal size of produced fluctuations is limited in 
this model] 2 ] even without back reaction. We have calculated the maximum value of (X 2 ) 
that would be achieved without back reaction for a wide range of values of q. It is plotted 
in Fig. by the dotted curve and for massless X has the label m x = 0. This curve can be 
approximated by (X 2 ) max oc q 12 , which is plotted in Fig. [T| by a solid straight line. 

Next, consider the case Mx ^ 0. The rapidly growing mass term of \ is a resonance 
suppressing factor. In the case of the A0 4 model, we will see that it works against the 
resonance extremely effectively. In the present model, however, the frequency of oscillations 
of the field ip is growing at the same rate as the effective mass of x, compare eqs. (§) and (||). 
In a sense, massive \ modes can remain tuned to the resonance, in some analogy with the 
conformal case. Still, for a non-zero Mx, the resonance terminates before q c s drops below 
unity. This is why the curves with m x ^ are shifted in Fig. [I] to the right with respect to 
the m x = case. 

In the absence of the Universe expansion and back reaction, the linearized equation for 
fluctuations, eq. would be a Mathieu equation. In that case, a necessary condition for 
the X particle production to fall into the strongest resonance band would be, for m x > 1: 
m x ^ g 1//4 , or equivalently q ^ (Mx/m) 4 . As one can see from Fig. [j], this estimate is four 
orders of magnitude off in the expanding Universe: e.g. with q ~ 10 4 the resonance barely 
develops even for m x = 1. To take the Universe expansion into account, we can replace q 
in the above estimate with g eff at some time before the end of the resonance stage, g e fr(' r i)- 
As evident from eq. (|Tl|), finding g e fr(Ti) is equivalent to finding 4>{ji). The quantities T\ 
and 4>(ti) in general depend on q. Although some analytical estimates of </>(ti) can be made 
||, its reliable determination requires our numerical results. We have found that for the 
range of q we consider in Fig. [I], one can use generically 4> 2 {ti) ~ 1CT 4 2 (O). Using q e s 
instead of q repairs the estimate: to efficiently produce particles of mass Mx, one needs 
<?efr(Ti) ~ (Mx/m) 4 , or q > 10 4 (Mx/m) 4 , in agreement with Fig. [l]. For Mx/m = 10, 
we get q ^ 10 8 , which for m = 10 13 GeV translates into g 2 > 4. x 1CT 3 . Considering the 
complexity of the phenomena involved, the estimate g 2 <; 10~ 4 made for this case in ref. (1^ 
is in a reasonable agreement with our result. 

The term "wide" parametric resonance is to some extent misleading. The resonance 
band itself is narrow, Ak g^ 4 . However, the resonance band can be anywhere in the 
interval < k TCS < g^ 4 , and indeed moves through this interval back and forth when q e g 
slightly changes. Now, since q e s is time dependent in expanding Universe, all k from this 



As opposed to the conformally invariant case P|T^R 
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range become resonant in turn|] and it is for that reason that one may consider q e g as a 
kind of effective width of the resonance. As a result of the motion of the resonance band, 
all k up to g 1//4 become filled up in the power spectrum, although to a varying extent. 

In Fig. [I], notice also the sensitivity of the maximal size of fluctuations with respect to 
small changes in q. Small changes in q can be simulated by small numerical uncertainties, 
such as that of the coefficient in the relation 0(0) = 0.28Mpi. It anyway makes sense to 
average out the rapid changes in Fig. [I] and concentrate on the general pattern. As a general 
trend, we observe a broad maximum of (A 2 ) max as a function of q, at fixed M X - It occurs in 
the region of q where the redshift and the back reaction play comparable roles in terminating 
the resonance. 

Now consider the limit when the maximum value attained by fluctuations is determined 
by the back reaction. This is the case of asymptotically large q. The general trend in 
(A 2 ) max in this region of q is given by 

<A 2 ) max ~ ^^M 2 , (12) 

which is plotted in Fig. [I] by the dashed straight line. In this region, we can also estimate 
(A 2 ) max in the Hartree approximation analytically, and the estimate is in rather good agree- 
ment with eq. fll2"|). It proceeds as follows. We first notice that because the decay of (f> 
into A is a stimulated process, it mostly takes place in a short period close to the end of 
the resonance, when A fluctuations are already large. Then, we should distinguish between 
the amplitude of the inflaton oscillations just before the decay, 4>(t\), which is determined 
mostly by the redshift, and its amplitude right after the decay, (frfo). For a wide resonance, 
ififa) can be smaller than 4>(t\) by orders of magnitude. Such a drop is clearly seen in Fig. 
|3] (even though for m x = 10 the value q = 10 8 is not an asymptotically large q). Because 
T2 — 7~i is such a short time, the Universe expansion during it can be neglected, and from 
energy conservation we get 

^(XVmW. (13) 

The second ingredient of our estimate is the observation [0 that the presence of (A 2 ) changes 
the effective mass of the inflaton; the effective mass is given by m 2 ff = m 2 + g 2 (X 2 ), where 
(A 2 ) is a slowly changing part of (A 2 ). When (A 2 ) becomes sufficiently large, m 2 ff w 
g 2 (X 2 ) ~ g 2 (X 2 ), and the effective resonance parameter becomes q c s = g 2 (j) 2 /Am 2 s ~ 
(p 2 1 (A 2 ). The resonance stops when q e s(r 2 ) ~ 1; using eq. (0), this gives[] 



It is perhaps not obvious a priori that this motion of the resonance band does not kill the 
resonance altogether. 

4 Our usual condition g e ff ^ (Mx/m e ff) 4 is satisfied for q e s ~ 1 at sufficiently large q, due to the 
large m e R. Notice also that because a large m e ^ lowers the Mx/m e B ratio, the Hartree correction to 
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The value of 0(ri) is weakly depends on q. As before, we use the generic estimate 
4>(ti) ~ 1O~ 2 0(O) for all q considered. We then obtain (X 2 ) max ~ 10 _3 Mp 1 / v /g, in rea- 
sonable agreement with the numerical result fljjD- A 1/ ' ^Jq (up to logarithms) estimate for 
(A 2 ) max appears in refs. |2]J^,P1|, although the way the estimate is obtained and the numer- 



ical coefficient are somewhat different from ours. We stress again that we do not know at 
present if the 1 / y/q behavior is an artifact of the Hartree approximation in this model. See 
also the discussion of the conformally invariant case below. 

Note that at asymptotically large q, the interaction term ^g 2 X 2 (p 2 dominates over the 
mass term of X during the entire resonance stage. So, at the resonance stage, the fluctuations 
of X cannot really be described as particles, since they are strongly coupled to the oscillating 
background. 

Model 2: Massless inflaton. In rescaled variables, the action for the matter fields in this 
model is 



S -I 
A 



(V^) 2 If. a \ 2 (V X ) 2 1 



2{ X -a X ) - — 



4 



Ml a 2 (r) 0^22' 



2A0 2 (O) a 2 (0) 



(15) 



where if and x are defined as before but r = v^A0(O)a(O)?7, £ = \^\(j)(0)a(0)x. 

We again take r = to be the first extremum of 4>(r). In this case we find 0(0) = 0.35Mpi. 
There is no resonance for a massive enough X (see below). In that case, after the inflaton 
oscillations start, the Universe quickly becomes radiation dominated and we can use 



o(r) 
o(0) 




0.51r+l (16) 



All other formulas are obtained from the corresponding ones for Model 1 by replacing m 2 
with \<fi 2 (0). Accordingly, we now use m x = Mx/\/A0(O). 

For a massless field x t ne problem becomes conformally equivalent to the one in the 
Minkowski space-time: the expansion of the Universe decouples in the conformal coordinates. 
The parametric resonance is always effective. We have done a complete non-linear analysis of 
this problem in the range 1 < q < 1000, using the approach developed in ref. 0. At q ^> 1, 
the evolution of the system is fast and we were able to trace the semiclassical thermalization 
0] up to a slowly evolving state with a developed Kolmogorov turbulence. These results 
will be described in more detail elsewhere; for previous discussions of kinetic phenomena 
after preheating, see refs. [||,|nj. For the limited set of values of q for which we now have 
the data, (X 2 ) at the end of the resonance stage in the full problem fits a 1/q curve. If 



the mass of <fi actually helps the resonance for a while, as we have indeed observed in our numerical 
results. 



S 



confirmed by more extensive calculations, such behavior would be in disagreement with the 
approximately 1 / y/q dependence obtained in the Hartree approximation. 

With a massive x fi e ld we have not found any parametric resonance for m x > 1 and q 
up to q = 10 6 . Analytically, the condition for the resonance to exist at time r (neglecting 
back reaction) is g 1 / 4 <: m x a(r)/a(0). In the conformally invariant case (m x = 0), the 
fluctuations reach appreciable strength at r ~ 100. Using that value in the estimate for 
m x 7^ 0, we would expect that for m x = 1, the resonance exists when q > 10 7 . We indeed 
observed a slow (physically insignificant) resonance at m x = 1 and q = 10 8 . 

In conclusion, we have studied the decay of massive and massless inflatons into massive 
excitations in the expanding Universe. We have found that a wide enough resonance can 
survive in the expanding Universe, though the expansion is very important for determining 
precisely how wide it should be. For a massive inflaton, the effective production of particles 
with mass ten times that of the inflaton requires very large values of the resonance parameter 
q, q ^ 10 8 . For these large q, the maximal size of produced fluctuations is significantly 
suppressed by the back reaction, but at least within the Hartree approximation they are 
still not negligible. We do not know at present if this conclusion will persist in the fully 
non-linear picture; the full account for rescattering is thus an important problem for the 
future. For the massless inflaton with a A0 4 /4 potential, the Universe expansion completely 
prevents a resonance production of particles with masses larger than y/\(j)(0) for q up to 
q = 10 6 . 
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FIGURES 




FIG. 1. Maximal value of the variance of fluctuations in Model 1 as a function of the resonance 
parameter. Dotted lines correspond to no back reaction and are fitted by a power law. Solid curves 
are obtained in the Hartree approximation. The dashed straight line is a power law fit to the 
Hartree data. 




FIG. 2. Fluctuation variance in Model 1, as a function of time, for two values of the mass of 
the fluctuating field; back reaction of fluctuations on the inflaton field was included in the Hartree 
approximation. 
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FIG. 3. Time evolution of the homogeneous inflaton field and the fluctuation variance in Model 
1 in the Hartree approximation, for given values of the parameters. 
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